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We study the effect of interaction on the temperature change in the process of adiabatic mixing 
two components of fermi gases by the real-space Bogoliubov-de Gennes (BdG) method. We find that 
in the process of adiabatic mixing, the competition of the adiabatic expansion and the attractive 
interaction make it possible to cool or heal the system depending on strength of interactions and 
the initial temperature. The change of temperature in a bulk system and a trapped system have 
been investigated respectively. 



Exciting development in ultracold atom systems has 
opened a new possibility to stimulate the many-body 
Hamiltonian that have been used to study strongly cor- 
related systems in condensed matter physics. Among 
the most exciting breakthroughs are experimental real- 
ization of the quantum phase transition from the super- 
fluid to Mott-insulating phase in bosonic system [l| , and 
the metal- insulator transition with fermionic atoms [2|, y] 
in optical lattice. Many of these many-body phenomena 
in ultracold atom systems are sensitive to temperature, 
for an example, a change in temperature due to adia- 
batic or non-adiabatic tuning of the parameters in the 
many-body Hamiltonian may hide the signature of the 
quantum phase transition, replacing it with a thermal 
transition instead. Adiabatic process, which keeps the 
entropy of the system a constant, is known to play an 
important role in experimental manipulation of the ul- 
tracold atoms, especially in cooling the many-body sys- 
tems. Therefore, how temperature changes in the pro- 
cess of adiabatic tuning the many-body Hamiltonian is 
a question of great interest and has important practical 
application in the experiment due to its potential rela- 
tion with cooling the ultracold atoms. In this process, 
the interaction have been known to play a key role in 
determining the temperature changing |4l4l3l|. and may 
lead to anomalous phenomenons 14 1. 



In this paper, we study the effect of interaction on the 
temperature change in the process of adiabatic mixing 
two components of fermi gases. Initially, the spinf and 
spinj, fermions are well separated, and the gases are com- 
pletely mixed finally. We introduce the s-wave interac- 
tion between the spin up and down fermions, which could 
be adjusted by s-wave Feschbach resonance, 15, 16]. We 
assume the interaction is attractive thus we can safely use 
real-space BdG method to analyze this question. With- 
out the interaction, the adiabatic mixing of the two com- 
ponent fermions would definitely cool the system. Be- 
low we focus on the effect of the interaction. Whether 
the process of the adiabatic mixing in this interacting 
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many-body system would cool or heating the system? 
We would show below that the answer to this question 
not only depends on the strength of the interaction, but 
also on the initial temperature of our system. 

The paper is organized as follows: firstly, we study 
the thermodynamic properties of an ideal situation and 
calculate the change of temperature from an ideal ini- 
tial configurations (a): a completely separate fermi gases 
to an ideal final configuration (b): a mixture of at- 
tractive fermions (as shown in FigH]). Further more, 
we propose an experimental setup to realize this adia- 
batic mixing process by introducing a gradient magnetic 
field, which recently been applied as a super cool atom 
thermometer [17] . We show that by tuning the gradient of 
the magnetic field from a large value to zero (as shown in 
FiglS]) slowly enough. Configuration (a) and (b) in FigH] 
can be connected adiabatically. Using the real space 
self-consistent Bogoliubov-de Gennes (BdG) method, we 
study the temperature change following the isoentrope 
and observe the variation of the real-space distribution 
of the particle number as well as the superconductor or- 
der parameters in the process of adiabatic mixing. In 
the final part of the paper, we briefly discuss the possible 
application and relations of our result on recent experi- 
ments about the mixing of two-component fermi gases. 

Firstly, we study an ideal situation. Considering two 
situations: one is the two-component fermions separated 
completely and there is no interaction between them, as 
shown in Figllja); The other is fermi gases fully mixed 
and the attractive interaction between spin up and spin 
down fermions lead to pairing between them, as shown 
in FiglTJb). Suppose there is a isoentropic process from 
Configuration (a) to Configuration (b), we address the 
question how do the temperature change? In the isoen- 
tropic process, it is well known that the expansion of the 
fermions would cool the system[lj,[l^, however, the mix- 
ing of attractive spin up and spin down fermions leads to 
pairing between them, which would heat the system con- 
sidering the entropy of the fermionic superfluidity is lower 
than the normal state in the same temperature. The tem- 
perature change obviously depends on the strength of the 
interaction since we are dealing with a many-body sys- 
tems. Further more, we find that it also depends on the 
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FIG. 1: Two configurations of two component fermi gases in 
liomogenous bulk system, (a) Separate noninteracting fermi 
gases; (b) Fully- mixed fermonic superfluidity. 



initial temperature of the well-separated noninteracting 
state. In the region of very low temperature, the pairing 
effect on the temperature change overwhelm that of the 
adiabatic expansion thus the temperature of the system 
would increase in this isoentropic process. When the tem- 
perature is close to the critical temperature Tc^ the pair- 
ing effect is not important and the temperature change 
mainly determined by the adiabatic mixing, which leads 
to the cooling of the system. 

For situation (a), the entropy and particle number at 
finite temperature in the non interacting fermionic gases 
are: are as follows: 



N^N^ + N^ = ^2kBl{e^-lXa,T) (1) 
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Sa = -2fc_B^[/(ek -Ma,r)ln/(ek - ^J.a),T) 

k 

+ /(-(ek-Ma),r)ln/(-(ek-/ia),T)] (2) 

while for situation (b), the corresponding thermody- 
namic properties are: 
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where /i^ is chemical potential and T is temperature, 
while S is the entropy of the system and N is the total 
particle number. f(x,T) is the function of fermi distri- 
bution at temperature T. a is the scattering length of at- 
tractive interaction. Ck = ^^^^-^k = \/(ek ^ A^b)^ + A^, 
where A is the pairing parameter for BCS state. The 
fermi momentums of free particle in two above situations 
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FIG. 2: (a) S-T curves in the homogenous bulk system 
for separated noninteracting fermi gases and fully mixed 
fermionic superfluidity. The interaction parameter is set as 
kjas = —0.6. Two dashed lines denote two isoentropic pro- 
cesses. /STi/Tp is the variation of temperature during the adi- 
abatic process (blue is cooling while red is heating process). 
(b)Dependence of the temperature change in the isoentropic 
process on the initial temperature T as well as the strength 
of the interaction. 



energy, respectively. The curves of entropy vs tempera- 
ture of such two situations are shown in Fig. 2(a). Ob- 
viously, the entropy of both situations must be extreme 
low near T = according to the third law of thermody- 
namics. 

We find fermonic superfluidity (BCS state) has the 
lower entropy while the temperature is not high enough. 
However, when we increase the temperature, the long- 
range order is destroyed gradually by the thermal fluc- 
tuation. When temperature has reaches to 0.75Tc, the 
entropy curves of attractive gases and that of completely 
separated gases intersect. If the temperature continues 
to rise, the fermi gases with attractive interaction change 
closer to normal mixture, of which the entropy is cer- 
tainly larger than that of the separated gas, as shown 
in our result. Cooper pairs in BCS state are completely 
destroyed at the critical temperature T,. of phase transi- 
tion. The curves of entropy has shown the competition 
between quantum order and thermal fluctuation. We as- 
sume the separated fermi gases can be mixed adiabati- 
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FIG. 3: (a) For the magnetic field with large gradient, 
fermions with different spin are well separated; (b) By adi- 
abatically tuning the gradient of the magnetic field to zero, 
the fermions are fully mixed. 



cally. The dashed lines in Fig. 2(a) denote two isoen- 
tropic processes with entropy fixed at 5*1 and ^2, and 
the first one is cooling while the seconde is heating. Fig. 
2(b) has shown the change of temperature in the adia- 
batic processes with different interaction strength, from 
which we can clearly see the heating and cooling regions 
respectively. 

Now we discuss the experimental realization of this 
adiabatic mixing process by applying a gradient mag- 
netic field, which was recently used as a super cool atom 
thermometer [171 . The gradient of the magnetic field is 
along the x-direction -B(x) = B{x)z. By tuning the gra- 
dient of the magnetic field from a large value to zero , 
configuration (a) and (b) in Fig|3]are connected adiabat- 
ically, as shown in FiglS] we assume that there is a hard 
constraint along x-direction while fermions are free along 
y and z directions : 
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To determine the thermodynamic properties such as the 
entropy and distribution of the particles in our system, 
we use the real-space BdG method to study this fermionic 
superfluidity under this gradient magnetic field. The 
hamiltonian in this case can be written as: 



H 



-h 



dxV&t (x)(^— V^ - /v + F.(x))vl>,(x) 



2m 



9 



dx*|(x)*^(x)*^(x)*t(x) 



(4) 



/icr is the chemical potential of spin a. Vo-(x) include 
trap potential Vtrap and Zeeman shift Vf'^'^. Here, we 
assume Zeeman shift of two components has form as 
V^'^^{x) = aXx (A > 0). Due to opposite energy shift, 
fermions with different spin are pulled to opposite direc- 
tion along X axis. Since there is not extra gradient po- 
tential and constraint along y,z directions, which means 
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FIG. 4: Distribution of density and pairing parameter. The 
total number of particle A'^ = 200 and the bare interaction 
parameter g is set as -1.937. We choose the temperature 
T = 2 X 10-^Tf, and (a)A=0, (b)A=0.016, (c)A=0.0286, 
(d)A=0.0446. 



the physical quantities in our systems such as density, 
fermionic pairing order parameter are uniform along y 
and z directions, thus they are only the function of x in 
real space. 

In order to show the density distribution in this in- 
homogenous system, we adopt to solve BdG equations 
in real space [18l - [20| . Within mean-field approxima- 
tion, the pairing gap and density are defined as A(a::) = 
— g(^j,(a;)^^(a;)), where g is the bare interaction parame- 
ter, ricix) = (^j^(a;)*cr(a;))- We use Bogohubov transfor- 



mation '^a{x) = ^i{uia{x)cio 
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the hamiltonian and get BdG equation: 
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and Ei-f = —Ei^, we can drop a in the above matrix 
equation. Pairing gap and density can be expressed 
by the quasi particle wave function ui{x) and vi{x): 
A(x) = <7Ei ui{x)vi{xYm), n^{x) = Y.i \ni{x)\' f{Ei) 
and n^{x) = I]j k/(^)l^ /("-^Oi where f{x) is Fermi- 
Dirac function at temperature T. We have adopt the hy- 
brid process introduced in [1^ [l^ . Besides the discrete 
spectra, the high energy part above an enough energy 
cutoff is included by local density approximation. The 
bare interaction g is replaced by a position-dependent 
effective interaction geff{x). Since the particle number 
of each component should conserve, we obtain the dis- 
tribution of pairing parameter A (a;) and density by self- 
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FIG. 5: (a) S-T curves in the trapped system with difTerent 
magnetic field gradient A = 0, 0.016, 0.0286. The bare interac- 
tion parameter g is set as -1.937. (b) contour of temperature 
change of the trapped system while the magnetic field gradi- 
ent is tuned from A = 0.0286 to A = adiabatically. 



consistent calculation. As shown in Fig. 4, where we can 
find fermions with different spin are separated gradually 
while the scale and amplitude of A (a;) decreases rapidly 
with the increasing of magnetic gradient. 



Now we return to address the same question pro- 
posed above: how do the temperature change in this 
adiabatic process. By using the equation of entropy 
S = -kBj:i{f{Ei)\nf{Ei) + f{-Ei)\nf{-Ei)], we can 
calculate the physical quantities and we can draw the 
curves of entropy vs temperature with different magnetic 
gradient A = 0,0.016,0.0286. (Notice that Fig.l(a) cor- 
responding two extreme condition of this case) Heating 
and cooling during the isoentropic processes are denoted 
by two dashed lines in the figuration. The change of tem- 
perature in the adiabatic mixing process as a function of 
initial state temperature T and strength of interaction 
g is shown in Fig. 5(b), where we set the gradient of 
the magnetic field is tuned from A = 0.0286 to A = 
adiabatically. As shown in the figure, both healing and 
cooling occur in our case. 

In conclusion, we have shown the change of tempera- 
ture for attractive fermi gases during the adiabatic mix- 
ing process in a homogenous bulk system and a trapped 
system. The competition between the effect of the adi- 
abatic mixing and the interaction on the entropy lead 
to interesting cooling and healing process. We focus on 
the attractive case where we can safely use mean-field 
method. Recently, the itinerant ferromagnetism has been 
discovered in a two-component fermi gases with repulsive 
interactions 12 1|. More recently, the mixing of two spin 
components of a strongly interacting Fermi gas have been 
realized experimentally [22i|. However, for the repulsive 
interaction, the mean-field result may be unreliable and 
the similar question how do the temperature change in 
the process of adiabatic mixing of fermi gases with re- 
pulsive interaction may be more interesting and deserve 
further investigation via other numerical method such as 
dynamic mean- field theory (DMFT)[23]. 
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